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Abstract 

We present a one-parameter family of exact solutions to Einstein’s 
equations that may be used to study the nature of the Green-Wald 
backreaction framework. Our explicit example is a family of Einstein- 
Rosen waves coupled to a massless scalar held. This solution may be 
reinterpreted as a generalized three-torus polarized Gowdy cosmology 
with scalar and gravitational waves. We use it to illustrate essen¬ 
tial properties of the Green-Wald approach. Among other things we 
show that within our model the Green-Wald framework uniquely de¬ 
termines backreaction for finite-size inhomogeneities on a predefined 
background. The results agree with those calculated in the Charach- 
Malin approach. In the vacuum limit, the Green-Wald, the Charach- 
Malin and the Isaacson methods imply identical backreaction, as ex¬ 
pected. 


1 Introduction 

In the Green-Wald framework [T] backreaction is represented by the effective 
energy-momentum tensor. It describes the leading-order effects of small-scale 
inhomogeneities on the global structure of spacetime. The main result of [T] 
states that this effective energy-momentum tensor satisfies the weak energy 
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condition and is traceless. Thus, the significant backreaction effect of small- 
scale inhomogeneities can be produced only by high-frequency gravitational 
waves and not by matter inhomogeneities. This result holds for all spacetimes 
satisfying the mathematical assumptions of Green and Wald [T] and, if ap¬ 
plied to cosmology, implies that backreaction cannot mimic the cosmological 
constant. The fundamental question is whether the Green-Wald framework 
captures the effect of inhomogeneities properly and whether it may be ap¬ 
plied to our Universe. Recently, this problem was the subject of a discussion 

0131 Eli- 

The most straightforward approach to clarify the issue described above is 
to search for exact spacetimes to which the Green-Wald framework may be 
applied and for which it leads to nontrivial backreaction effects. The study 
of such an example is the main aim of our paper. 

The theorems of Green and Wald [1| are based on mathematical assump¬ 
tions that are not intuitive. The main assumption is the existence of a 
one-parameter family of solutions to Einstein’s equations, g a b( A), A > 0, that 
approaches the effective (background) spacetime in a special way as A —* 0. 
It is not necessary to know this family in an exact form to predict the prop¬ 
erties of the backreaction effect. However, one should provide examples of 
such families to explore the class of solutions to which the Green-Wald frame¬ 
work applies. Moreover, since there exist other similar frameworks, one may 
use these examples of one-parameter families of exact solutions to check the 
consistency of different methods. In this paper, we compare the Green-Wald 
framework to the Charach-Malin approach [6] and, in the vacuum limit, to 
the Isaacson method m- Although these methods do not differ much when 
it conies to practical calculations, they are based on different mathematical 
assumptions. The Charach-Malin and Isaacson methods do not involve ‘tak¬ 
ing limits’, and thus there are no doubts that they describe backreaction for 
finite-size inhomogeneities0 The essential advantage of the Green-Wald for¬ 
malism over other similar approaches is the fact that it allows one to prove 
general theorems about the properties of backreaction. These theorems do 
not have their counterparts in other approaches. Therefore, it is interest¬ 
ing to study relations between the different frameworks. The exact solution 
presented in this paper provides such an opportunity. 

1 This aspect of the Green-Wald method has been recently analyzed in [5]. 


2 



Some time ago w^] presented a one-parameter family of exact solutions 
to Einstein’s equations that was used to test the Green-Wald backreaction 
framework [9|. This was the first nonvacuum example of a family of ex¬ 
act solutions that satisfies the Green-Wald assumptions. The simplicity of 
this family makes it an interesting testbed for this formalism. On the other 
hand, the energy-momentum tensor in |9j is uniformly convergent as A —» 0, 
so its weak limitj^] coincides with the ordinary limit—such matter inhomo¬ 
geneities do not model the inhomogeneities in our Universe sufficiently well£j] 
Therefore, it would be instructive to find a nonvacuum example such that 
the weak limit of the energy-momentum tensor is nontrivial and the density 
contrast of the matter is high. We provide such an example in this paper. In 
our example, the energy-momentum tensor and the Ricci scalar are not uni¬ 
formly convergent as A —» 0. The density contrast of matter does not vanish 
and remains bounded, but the frequency of the density fluctuations blows 
up as the background spacetime is approached. We show that in our model 
this kind of inhomogeneities in the distribution of a massless scalar field do 
not contribute in the leading order to the backreaction effect in accord with 
what was proved in a general setting in |:f|. The bounded nonvanishing am¬ 
plitude, high-frequency matter inhomogeneities studied in our article provide 
an example of an intermediate behavior between decreasing amplitude, high- 
frequency matter inhomogeneities studied in [S] and unbounded amplitude, 
high-frequency matter inhomogeneitie^J an example of which has yet to be 
found. 

There already exist four articles in which Green and Wald presented their 
method HI 0 US HI] and we refer the reader to those papers for mathemat¬ 
ical details. However, we think that it could be instructive to review their 
framework from a different perspective. 

2 In collaboration with Krzysztof Glod and Alicja Konieczny. 

3 For a definition of a weak limit see (T] . 

4 This was pointed out by Stephen Green and Robert Wald and suggested in the sum¬ 
mary of the paper coauthored by us [S]. Recently similar remarks were presented in [2] - 

5 This classification of inhomogeneities is based on the behavior of the amplitude and 
frequency of fluctuations of matter density as the background spacetime is approached 
along the one-parameter family of solutions. 
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2 Averaging and the Green-Wald framework 

In this section we briefly review the Green-Wald framework. We restrict 
ourselves to the ‘leading-order analysis.’ We keep the mathematical material 
to a minimum, however, for the sake of clarity it seems necessary to introduce 
some definitions. For brevity, sometimes the surnames Green and Wald will 
be abbreviated as GW. 

Let M be a Lorentzian manifold and L(M) a set of (at least) C 2 Lorentzian 
metrics on M. Following Green and Wald, one may introduc^] the no¬ 
tion of metrics with small-scale inhomogeneities. We denote this subset of 
L(M) by Li n h(M). Hence L e ff := L(M)\Li n h(M) corresponds to all metrics 
without small-scale inhomogeneities. Any element of this set, if not indi¬ 
cated otherwise, will be denoted with g and called an effective metric. 
By g\ C L in h(M) we define a one-parameter family of metrics, such that 
= lim A ^ 0 + ,9a £ L e //(M) (uniform convergence). 

It is usually assumed that the backreaction problem will be solved if an 
appropriate averaging procedure is found. This procedure—let us denote 
it by srf : L in h{M ) —> L e ff(M )—if applied to an inhomogeneous metric g, 
should give an averaged effective metric and thus &/(g) = fy°b Then, it 
is natural to expect that a well-defined stf will satisfy some basic properties: 
it should be covariant and unique. 

Unfortunately, it seems that this kind of approach is not the best way 
to proceed in most physically interesting situations. For example, finding a 
correct sf is not sufficient to solve the cosmological backreaction problem: 
one still needs ‘the true metric of the Universe Q— let us denote this hypo¬ 
thetical metric as gu G L mh {M^ — to average it to find the effective metric. 
It is usually expected that srf(gu) = gnw, where gmv G L e //(M) stands 
for the Robertson-Walker metric. However, if we had known gu, then most 
important problems would be solved and one would need the effective metric 
g RW only to simplify calculations. On the other hand, although knowing srf 
alone does not solve the problem, it could be helpful. One might apply srf to 
some toy models, study how small-scale inhomogeneities alter the averaged 
effective metric and try to generalize the results to a cosmological context. 

6 For example, as presented in the Introduction in jB] . 

'Or, maybe better, ‘the metric that encapsulates all essential properties of the true 
metric of the Universe.’ 

8 We assume here that M corresponds to a manifold which may be used to construct a 
model of the Universe. 
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The Green-Wald framework that we are going to study in this paper was 
not devised to work in the way described above. This approach is a gen¬ 
eralized perturbation theory and as in every perturbative approach, i.e. the 
Isaacson high-frequency limit [Hi, full covariance is lost. Since it could be 
not clear what a perturbative approach may mean in the context of backre- 
action, we clarify this issue below. 

First, as in every perturbative approach, the effective spacetime must 
be assumed a priori. This is not a disadvantage of the framework, because 
the theorems of Green and Wald [1] hold for any background spacetime and 
arbitrary ‘perturbations’ satisfying their mathematical assumptions. Thus, 
the properties of backreaction in this framework are known in general and a 
particular form of the background metric was not needed in the proofs. The 
GW method is in some sense an inverse approach to the one presented at 
the beginning of this section. If srf is an arbitrary averaging procedure such 
that &/(g) = g(°\ then C L(M ) corresponds to the whole class 

of possibly inhomogeneous metrics that may be effectively represented by 
g(°\ Similarly, in the GW approach ^gw{9^) corresponds to a set of one- 
parameter solutions to Einstein’s equations g\ such that lim / \^. 0 + 9x = g ^ 
(uniform convergence) and that satisfy the remaining assumptions of Green 
and Wald [(iii) and (iv) in [I]]. |Hereafter, the Green-Wald assumptions 
(i)-(iv) will be termed for brevity the GW assumptions.} Thus, the GW as¬ 
sumptions provide a precise definition of s^Q^{g^) which corresponds to a 
‘cloud’ of possibly inhomogeneous metrics that are in a wea/0 sense not far 
from some effective background g(°\ In the GW framework essential infor¬ 
mation about inhomogeneities is encoded in first and second derivatives of 
metrics in This idea is best known from the Isaacson approach^] 

How is backreaction defined in this kind of perturbative frameworks? The 
metrics in ^owig 1 ' 0 ' 1 ) and the corresponding energy-momentum tensors sat¬ 
isfy Einstein’s equations. We would like to represent these inhomogeneous 
metrics by an effective metric g^ and a smoothed energy-momentum ten¬ 
sor T^ 0) {^] However, there is no reason why both quantities should satisfy 
Einstein equations. We still may pretend that they do if we add an addi- 

9 The precise meaning of this word follows from the GW assumptions [Tj. 

10 The Isaacson approach was originally applicable only to vacuum spacetimes. It was re¬ 
formulated in a mathematically rigorous way by Burnett m- The Green-Wald framework 
(Tj generalizes the Burnett result to nonvacuum spacetimes. Nonvacuum direct general¬ 
izations of the Isaacson formalism are also possible m- 
11 Green and Wald provided a procedure to determine . 
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tional term, so G(g^) = 87 r(T^ 0 ^ + f(°)). This term is called the effective 
energy-momentum tensor. It represents the effect of inhomogeneities on the 
background spacetime. So, as it possibly happens in cosmology, if one would 
naively construct a cosmological model with g(°\ then one would discover 
that the corresponding energy-momentum tensor contains a strange contri¬ 
bution U°h This problem is known as a ‘fitting problem’ n. One assumes a 
priori that Einstein’s equations hold for the effective metric g^j^jand then 
asks whether the corresponding energy-momentum tensor contains some un¬ 
expected corrections in the form of We point out that this kind of inverse 
approach is usually adopted in modern cosmology, where philosophical prej¬ 
udices supported by astronomical data were used to construct the standard 
cosmological model. With the new observations the ‘best fit’ model is being 
improved, but one roughly knows the effective cosmological spacetime. The 
question is whether the parameters that define it have been interpreted cor¬ 
rectly. In particular, is the energy-matter content what it seems or is a part 
of it an artifact of averaging? 

We have mentioned at the beginning of this section that one should not 
see s^gw as a straightforward averaging procedure. Indeed, if s^gw is applied 
in this way, it acts on one-parameter families of metrics g\ which satisfy the 
GW assumptions, but its action is trivial: s^gw — b m A-s>o+- It seems to be 
more appropriate to define the GW framework as a map : g\ —>■ where 

g\ satisfies the GW assumptions, and thus already encodes information about 
the effective metric = liniy_ >0 + 9\- This allows one to prove some general 
theorems about relations between inhomogeneous metrics £^Gw(g^ 0 ' 1 ) and the 
energy content of the effective spacetime with the averaged metric g^ which 
represents all these inhomogeneous metrics [Tj. Since there was no need to 
assume a particular form of the effective metric g(°\ Green and Wald’s results 
tell us about properties of backreaction in general. Using notation introduced 
above the main result of Green and Wald |1] may be formulated as follows. 

Theorem, (by Green and Wald JVj, 2011) 

Let g\ be a one-parameter family of solutions to Einstein equations satisfying 
the GW assumptions. The effective energy-momentum tensor t^ := ^(gx) 
satisfies the weak energy condition and is traceless. 

This theorem implies that only gravitational radiation may contribute in 
the leading order to the backreaction effect, and thus in cosmology small- 
scale inhomogeneities cannot mimic a cosmological constant. This theorem 

12 In cosmology, this choice is based on a best fit to observational data. 
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also shows the advantage of the Green-Wald approach over the Isaacson-like 
nonvacuum extension of the high-frequency limit |6lfl5] and justifies the need 
to introduce one-parameter families of metrics g\. In the original Isaacson 
approach (vacuum spacetimes) positivity of the effective energy density was 
shown only under an additional ansatz The reason for that is a limited 
control over the coordinate dependence. Thus, the GW approach is not only 
a mathematically rigorous extension of the Isaacson approach, but there is a 
qualitative difference: within the GW framework one may study the nature 
of backreaction itself (not only backreaction in particular models). On the 
other hand, this kind of approach also has some disadvantages. Imagine that 
we know gu and we expect that it may be effectively represented by a g R w 
metric. It would not be easy to estimate the size (not nature, but size) of 
the backreaction effect within the GW approach^] One would need to use 
observational data to determine the appropriate member of a g R \y class and 
find a family of exact g\ metrics such that g\ = gu, lim A _>. 0 g A = g RW In 
addition to that, there may exist in principle many such familie^] g\ so in 
such a case, for a given g R w and gu, the effective energy-momentum tensor 
t would not be uniquely specified]^] Thus, the power of the GW approach 
comes from the rigorous identification of a ‘cloud’ of inhomogeneous metrics 
that are in some sense not far from the effective metric. This makes the 
general studies of the nature of backreaction possible. The price for this 
identification has the form of the GW assumptions about the existence of 
appropriate one-parameter families which limits the ability of the framework 
to provide quantitative predictions. We are not aware of any other solution 
of this identification problem which would give easy and precise control over 
coordinate dependence] 17 ] 

Finally, we return to the fundamental question of whether the Green- 
Wald framework applies to our Universe. It seems not to be easy to identify 
elements of ^Qy^(gRw) and it is not obvious if the ‘true metric of the Uni- 


13 Since gravitational waves have been detected m , it is of great importance to estimate 
their cosmological backreaction effect. Although we have good reasons to assume that it 
is negligible, one may still insist on calculating its precise magnitude and this aim is hard 
to achieve within the GW framework alone. 

14 The parameter A may always be rescaled to obtain gu for A = 1. 

15 This is called a ‘path dependence’ in [2.. 

16 We study this issue for a particular one-parameter family g\ in Subsection 

17 In the GW approach any one-parameter family of diffeomorphisms 'Iq is allowed as a 
coordinate transformation provided that 'f'o = id (or more generally: any 'U that modifies 
the limit lim^o g\ trivially is allowed). 
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verse’ gu belongs to this class. 18 The most straightforward approach that 


may clarify this issue is to look for any exact families of metrics which be¬ 
long to £^Gw(g < ' 0 ' 1 ) for some effective metric g which is not necessarily of 
physical interest. A priori there is no reason for the GW assumptions to be 
very restrictive^ however, there are several logical possibilities to explore]^] 
Any of the conditions described below would invalidate the GW framework as 
a proper tool to describe the backreaction effect of matter inhomogeneities: 

1. ^civ(fi ,(0) ) = 0 for an y 9 {0) - 

It follows from examples presented in (9j 10] [12] that this is not the 
case. 


2 . &? G w contains only vacuum metrics for any g(°\ 

This possibility was excluded by the example presented in [9]. 

3. £^Gw(g ( ' 0 ' 1 ) contains only metrics with trivial matter inhomogeneities 
and gravitational waves on the g ^ background. 

Such a possibility would naturally explain why the effective energy- 
momentum tensor is traceless. The example presented in this paper 
excludes this scenario. 


4. gu e ^gw(9Rw) an d 3 g\ satisfying the GW assumptions such that 
9i — 9Ui lim^o , 9 a = Qrw for the observational best-fit gnw, but &(g\) 
does not give the correct 


The possibilities 1 and 2 seem unnatural and they have already been ex¬ 
cluded. The aim of this paper is to exclude the possibility 3. In this article 
we present a one-parameter family of solutions to Einstein’s equations g\ 
which satisfies the GW assumptions and which has nontrivial matter inho¬ 
mogeneities. By the ‘nontrivial matter inhomogeneities’ we understand in¬ 
homogeneities which are represented by the energy-momentum tensor which 
is not uniformly convergent in the limit A —> 0. 

18 One may construct toy models that in some sense exhibit the backreaction effects, but 
cannot be studied directly within the GW framework mm- 

19 What is really restrictive is the fact that we are looking for solutions that are exact 
and, in addition, have a concise mathematical form. 

20 Since h°) is traceless one should check if any solutions with nontrivial matter inhomo¬ 
geneities may belong to = 0 ? gw(9 ^)- Otherwise, h°) being traceless would be an artifact 
of assumptions that are too restrictive. 
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3 Einstein-Rosen waves and a massless 
minimally coupled scalar field 

The cylindrically symmetric metric with two hypersurface orthogonal Killing 
fields dz, dp may be written in the form 

g = e 2 ^~^ (- dt 2 + dp 2 ) + p 2 e~ 2i 'dip 2 + e**dz 2 , (1) 

where p > 0 , — oo < t,z < oo, 0 < p < 2 tt and the metric functions ^ 
and 7 depend on t and p only. Some nontrivial vacuum spacetimes with 
the metric (0 were discovered by Beck ESI, but they are better known as 
Einstein-Rosen waves [201 EH- In this paper, we will investigate nonvacuum 
generalization of these solutions—a massless minimally coupled scalar field 
0 (f, p) will be added. 

There is an almost direct correspondence between generalized Einstein- 
Rosen waves and generalized Gowdy cosmologies. A subset of solutions stud¬ 
ied in this paper may be reinterpreted as polarized three-torus Gowdy models. 
Thus, the example that we provide may be seen as a nonvacuum generaliza¬ 
tion of the Gowdy solution presented by Green and Wald in [10] and a special 
case of the polarized Gowdy cosmologies investigated by Charach and Ma- 
lin [6]. Since Charach and Malin proposed their own approach to study the 
high-frequency limit of their solutions [6], we find it instructive to compare 
our results. This is done in Section E 

The energy-momentum tensor of a massless scalar field has the form 

Tab = d a (f)d b (t) - ^g a bd c (t>d c 4) . 


It is convenient to rescale the scalar field </>, so from now on we will use 0, 


where (f> = Einstein’s equations reduce tcpj 


ip" + V - ij = 0 , 

P 

(2) 

y = p (<y+y 2 +y 2 +y 2 ), 

(3) 

7 = 2 p . 

(4) 

21 In the following, dots and primes denote derivatives with respect to t and p, respec¬ 
tively. 
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The scalar field 0 satisfies 


0" + V - 0 = 0 (5) 

P 

which corresponds to V a V a 0 = 0 and follows from the contracted Bianchi 
identities. The equations ([2]) and (|5]) are linear and identical to the equations 
for a polarized cylindrical wave propagating with the speed of light in Eu¬ 
clidean space. These equations are decoupled and may be solved separately. 
Once solutions to the equations (|2]) and (J5| are known, the remaining metric 
function 7 may be found via the equations (|3]) and Q using quadratures. 
Moreover, the form of equations (|3]) and Q implies that every smooth pair 
of functions 0, 0 give rise to a smooth metric ([Tj) . The Minkowski space- 
time corresponds to 0 and 0 being constant. Decoupling of the dynamical 
equations for 0 and 0 implies that every solution of the vacuum equations 
(0 = 0) may be trivially extended to contain a nontrivial scalar field. [For 
a given 0, it is sufficient to take any 0 satisfying (|5| and calculate the new 
function 7 via quadratures.] 

In the succeeding section we will need the energy density of the scalar 
field e as measured by the observers comoving with the coordinate system 
(with the four-velocity u = e^~ 1 d t ) 

e = T ab u a u b = ^e 2 ^ (> + 0' 2 ) , (6) 

and the Ricci scalar R, which, with the help of the equations ([2]) and (|5]) , 
may be written in the form 

R = 2e 2W ’" 7) (V 2 - 0 2 ) • (7) 

In general, spacetimes considered in this section may contain gravitational 
radiation. The analysis of backreaction of high-frequency gravitational waves 
in these spacetimes was presented by Podolsky and Svitek [13]. 

4 One-parameter family of solutions 

Our construction of a one-parameter family of metrics satisfying the assump¬ 
tions of the Green-Wald framework is based on the solution described in the 
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previous section. Namely, we choose the following particular solutions of the 
equations ^ and Q, respect ivelyp^j 

(f>x{t,p) = aVX F x (t,p) , i>\(t,p) = {3\/\ F x {t,p) , A > 0 , (8) 


where F\(t, p) = Jo (|) sin (|) and A is a parameter. J 0 is the Bessel function 
of the first kind and zero order. The constants a, f3 are real and indepen¬ 
dent of A. Integrating the equations ([3]) and Q and setting the additive 
integration constant to zero we obtain 


7a (T p) 


(« 2 + n 2 

2A P 


“'» ( A ) + J?( A 


Ps 



Our one-parameter family of solutions to Einstein’s equations, denoted g a b{ A) 
(where A > 0 is the parameter), has the form ([TJ) with the metric functions 
t/b 7 and the scalar held given by ([8]) and 0- 

For p/X » 1 the asymptotic behavior of the Bessel functions is given by 


<> 



cos(p/A 



+0 $ 


( 10 ) 


Therefore, in the limit A —» 0, we have 

-t 0 , 7a (cn 2 + /3 2 )p/ir , </>A -t 0 . (11) 

The background metric := limy-^ gab(X) is curved and has the form 

g(0) = e 2{a*+p) P /n ^_ dt 2 + dp 2 j + p 2^2 + d ~2 _ (12) 


The functions ( [IT| do not satisfy the equation (j3|, and hence the metric g (0 ^ 
does not belong to the class of solutions described in Section [3j 

Let A\(t, p ) = Jo (y) cos (|) and B\(t,p ) = Ji(^)sin(|). Then for 
A > 0 the nonzero components of the energy-momentum tensor of the scalar 

22 0ne may consider a more general class of solutions, but we would like to make our 
example as simple as possible. 
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field (f>\ take the form 


rr 2 

Ttt( A) = T pp (X) = —(Al + Bl), (13) 

Tt P ( A) = T pt (\) =-~A X B X , 

r w ( A) = ^x e ~ 2lx P 2 ( A l- B D’ 

T- z - z { A) = T w (A)p- 2 e 4 ^^ . 

The energy density of the scalar held measured by a comoving observer ([6]) 
is 

€(A) = + Si) ■ (14) 

The Ricci scalar (Jr]) equals 

R( A) = 2(L e 2( /5VA^-7A ) ( jB 2 _ ^ _ (i5) 

A 


Since the Bessel functions are regular at zero, it follows that the energy 
density and the Ricci scalar remain bounded as p —> 0 for any solution with 
nonzero A. 


In the formulas (13), (14), and (15) the auxiliary functions A \, B\ may 


be approximated for a small value of A (p/A 1) with the help of (10) by 


A x (t,p) « 

2X 

\ — cos 
V vrp 

a- 

^ 1 ^ 

1 

) cos(^) , 

(16) 

B\ (t,p) « 

\ll- p sin{ 

f p 
vA 

4/ 

1 sin(^) . 

(U) 


In this approximation, the singular factor 1/A in T a b, e, R cancels and these 
quantities exhibit rapid oscillatory behavior with a finite amplitude of os¬ 
cillation^]— the case which we refer to as bounded nonvanishing amplitude, 
high-frequency inhomogeneities. 

23 The frequency of oscillations blows up as A —> 0, but the amplitude of oscillations is 
bounded and nonvanishing in this limit. 
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5 Inhomogeneity effect 


The equation satisfied by the background metric has the form [Tj 

G ab (g {0) ) = 8tt T® + 8tt , 

_ofT ab (A) 

T ab = w-limT a6 (A) • 


where Ty 1 is so-called the weak limit 24 


A-s>0 


(18) 


The new term on the right-hand side, denoted as t a ° b \ is called the effec¬ 
tive energy-momentum tensor and it encodes the backreaction effect of in¬ 
homogeneities. This tensor represents the additional terms that arise in 
the averaging process of Einstein’s equations, so that these equations may 
still hold for A = 0 with the modified energy content. A straightforward 
way to determine the effective energy-momentum tensor is to calculate 


(o) 

ab ' 


€ = £<W») - T c 

It follows from (13), (16), and (17) that T a b( A) is not uniformly convergent 
as A —> 0. This is an essential novel property in comparison to the example 
published in [9j. In the weak limit the nonzero components of the scalar field 
energy-momentum tensor are 


rp(®) _ 'WO) _ _ 

~ pp ~ 8ir 2 p ’ 

so in the investigated model a fast-varying scalar field may be approximated 
by an anisotropic null fluid. Since the nonzero components of G a b(g ( '°' > ) are 


G it (g <0) ) = G„(g (0 >) 


a 2 + /? 2 

7rp 


then the nonzero components of the effective energy-momentum tensor are 



(19) 


The effective energy-momentum tensor above is traceless and satisfies the 
weak energy condition, as predicted by theorems in [1]. It also has the 

24 We use the symbol w-lim to denote the weak limit. 
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form of an anisotropic null fluicj^] and depends on (3, but does not depend 
on a (the parameter a controls the magnitude of the scalar field). Hence, 
inhomogeneities of the scalar field do not contribute in the leading order to 
the backreaction effect. 

The Green-Wald framework requires the existence of a smooth tensor 

£elc0 


Habcdef ■ W lilU [V a /r cc ;(A)V^/? e j(A)] , 


where h ab ( A) 
be written in 


8 ^ 2 ’ = 


:= 9ab{ A) — f/a b ■ The effective energy-momentum tensor may 
terms of Pabcdef 



n c d 
r 1 ca bd 


( 20 ) 


Although we have already calculated from ( fl8| , it is instructive to de¬ 
termine habcdef and check the consistency of the GW framework. A lengthy 
calculation performed with the help of MATHEMATICS^] yields 


fttttttt ftttpppp ftttttpp 

ftpptttt ftpppppp ftppttpp 

9 

" n2 3 

fttt<p<p<p<p ft ppip^pip^p ~PP 5 


-fr't A(« 2 + /3 2 ) : 

nr x 7 


7r 


3 4(a 2 +/3 2 )p/7r 


ftttz 


— — Z. R 2 _1 

ftppzzzz P p 1 


7r 


ftttpp(f<p ftttttcpip ftppppipip ftpptt(f(f ~/3 /9 6 ^ ^ ? 

2 2 2 

ftttppzz ftttttzz ftpppp zz ftppttzz — )p/ \ 

7T 

__2 2 

Httififzz Hppipipzz P P ■ 

7T 


25 Similarly an energy-momentum tensor containing terms of the form of an anisotropic 
fluid appears in another approach to averaging [25]. 

26 The derivative operator V a is associated with the background metric g^. 

27 Some calculations were carried out with the xAct package for Mathematica |23j . 
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All other components follow from symmetries [T] 


M abode f l^(ab)(cd)(ef) H'abefcd 


or are equal zero. We substitute this tensor into (20) and calculate compo- 
nents of t^j. They coincide with (19) in accord with the general calculations 

in pq. 


6 Generalized three-torus Gowdy cosmologies 


The solutions studied in this paper may be reinterpreted as generalized 
Gowdy cosmologies. The main aim of this section is to compare our results 
to the high-frequency limit presented in [6|. 

Charach and Malm found a general solution which represents polarized 
Gowdy cosmologies with three-torus topology that are minimally coupled to 
a massless scalar held [6]. Using their notation the metric may be put into 
the following form: 

g = L 2 [e^(—d£ 2 + dz 2 ) + ^e p dx 2 + ^e~ p dy 2 ] , (21) 


where L is a constant (for simplicity we assume from now on that L — 1), 
0 < z < 2 tc, £ > 0, 0 < x,y < 2 tt. The metric functions /, p and the scalar 
fielcp^ ip depend only on z, £; periodicity in z is assumed. The following 
formal complex substitution of variables and redefinition of metric functions 
brings the metric Q into the form: (21) 




(■ iz,i£,iy,x) , 

Tp(t,p) 

-> 

+ !n£) , 

7 (t,p) 

—>• 

+P(z,0 + ln£) , 


-> 

2^r(p(z,£) . 


( 22 ) 


The variables z and £ are assumed to be real, and thus one generates a new 
solution to Einstein’s equations from the old one by a complex substitution. 

28 The letter ip denotes the angular variable in our original metric 0 (not a scalar 
field), but we prefer to keep the original notation as in [5] . Hopefully, this will not lead to 
confusion. 
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We note that the original metric ([!]) was not periodic in t, but the metric (21) 
is assumed to be periodic in 0 with the period 2i r. In addition, it is assumed 
that x is a periodic variable, but z was not periodic for Einstein-Rosen waves. 
The original range of coordinates t, z has changed. In three-torus Gowdy 
models they correspond to z and x, respectively, and are periodic with the 
period 2n. This implies that the cosmological solutions must be periodic in z, 
so the functions ^(t, p), 7 (t,p), 4>(t, p) are substituted by periodic functions 

p( z ,0 = p( z + 27rfci,f), f( z ,0 = f( z + 27rfc 2 ,0, <P( Z >0 = P( z + 27rfc 3 ,0> 

where ki are arbitrary integers. We also note that time and space coordinates 
are swapped: p was a spatial coordinate, but £ is a time coordinate and t 
was a time coordinate, but z is a spatial coordinate. 

The solutions ^ and 0 studied in this paper are related to a general 
solution discovered by Charach and Malin [ 6 ]. Since Charach and Malin 
assumed periodicity in z (three-torus cosmologies), and we did not assume a 
periodicity in t, this correspondence is not one to one. Let j be a fixed large 
integer j>l and choose A = 1 /j. (In some calculations A will be treated as 
a discrete parameter that approaches zero as j —» + 00 .) In order to see the 
correspondence, one should substitute (no summation over j below) 


«0 = -1 , An = » C n = ^-/=^n a /V~j , * 

2v 7r 


= — s j 

2 j n 


(23) 


into solutions (16), (17), (18), and (19) in | 6 ], set j = 1/A and put all the 


remaining constants to zero. Finally, the inverse transformation to (22) gives 


and (J9]) from (16), (17), (18), and (19) in [6]. [Note that the formula for 
the function fs was not described properly in [Bj. In order to obtain it from 
their equation (19), it is necessary to replace ao by /3o, multiply all terms by 
a factor 167 T and replace A n by CJ n and B n by D n .\ 

Using the procedure above we find that our solution after the substitution 


( 22 ) corresponds to the metric ( 21 ) with 


p = -Inf+ 2/3\/AJi,(|)sm(^) , 


^ = 27F“' /X/ 4 )sill( V 


a 2 + /3 2 


A 


/ = i - £ 2 J o(\) + 


'A y 




'i' 

'A y 


.z, 

y 


-2^VA/ 0 (|)sin(p , 


( 24 ) 


where ip denotes the scalar held. 
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6.1 The high-frequency limit 


It is instructive to calculate the high-frequency limit of our solution using 
the Charach-Malin procedure and compare it to the high-frequency limit in 
the Green-Wald approach. 

Using (10) for £/A 1 we obtain to leading order in A/£ 


2a/2 /?A /£ 7tA . z 

p “ _ + t? y? cos (a _ ij sm( A ) ’ (25) 

1 aX /£ tt\ . ,z s 

* " 777f COS lj“4 / ) Sm( A ) ’ 

/ ss !(a 2 + / 3 2 )t, 

7r 


in accord with the formulas (34), (35) and (36) in |6] (with our choice of 
constants). Following Charach and Malin we decompose the metric (21) into 
the ‘background’ g and ‘wave’ part h as g ~ rj + h, where 


T] = e 2(« 2 +/3 2 )€A(_^ 2 + dz 2 ) + ^2 + , ( 26 ) 

h = p(dx 2 — £ 2 dy 2 ) , 
p = p + In £ . 


We assumed that p is small. After appropriate redefinition of variables 
and functions Charach and Malin’s background metric rj corresponds to the 
Green-Wald background metric g ^ given by (12). 


In order to show that the background geometry is created partly by the 
gravitational waves and partly by the scalar held Charach and Malin cal¬ 
culated the energy-momentum tensor of the background metric r/ and de¬ 
composed it into the traceless part T 1 - 1 ) and nontraceless part T hence 
T = jA 1 ) + JA 2 ) = g Y -G(rj). For our solution we find 


t = T m = + <fe 2 ). r (2) = o, 


(27) 


where the nontraceless part T vanishes because of our choice flo = 0 (/3 q is 
one of the constants in the full Charach-Malin solution). The traceless part 
T W is further decomposecQ into the gravitational waves part T GW and the 

29 This decomposition is based on the interpretation of constants, e.g. for our solution 
a controls the amplitude of the scalar field y>. 
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scalar radiation part T sn 


^( 1 ) _ r-pGW rp 


SW 


jiGW 


/ 3 2 


a 


8tt 2 ^ 


(de + dz 2 ), T = + dz 2 ) ■ 


Let be a null vector defined by its covariant components as 

k = _ (d^ + dx) . 

V *£ 

Then, 8nT^ = f3 2 k^k u and 8ttT^ v = a 2 k ll k u . Thus, Charach and Malin 
interpreted both of those terms in the traceless part of the energy-momentum 
tensor as a null fluid representation of collisionless flows of ‘gravitons’ and 
scalar massless particles. This is one of the main results in their paper [6]. To 
sum up, we started with an inhomogeneous cosmological model with small- 
scale matter and ‘gravitational field’ inhomogeneities and showed how to 
represent this model by an effective spatially homogeneous (but anisotropic) 
spacetime filled with collisionless flows of massless scalar particles and ‘gravi¬ 
tons.’ 

Up to this point calculations in the Charach-Malin framework were al¬ 
most identical to those in the Green-Wald approach, so it is a good moment 
to summarize the differences. First, Charach and Malin split the metric into 
rj and h using their intuition and the interpretation of the metric r). In fact, 
this is done via ‘inspection’ of the line element. In the coordinate system that 
they have chosen, the form of the line element has a natural split into the 
‘background’ and the ‘wave’ part. This gauge-dependent intuition is strongly 
supported by a nice physical interpretation of T — d-G(r)); a priori there is 
no reason why T calculated in such a way should correspond to any physi¬ 
cally reasonable form of the energy-momentum tensor. In the Green-Wald 
approach is calculated as the limit A —>■ 0 of the one-parameter family 
of metrics g( A). Similarly to the Charach-Malin approach, we must define 
in advance our background metric. Both backgrounds coincide because the 
coefficients A n and C n that we have chosen approach zero in an appropriate 
way in the limit A — > 0. We are interested in the backreaction effect of finite- 
size inhomogeneities (for a finite 0 < Ao 1) and up to this point there 
would be nothing ‘ultra-local’ in our calculations if we had conducted them 
in the Green-Wald framework. The result would not depend so far on the 
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choice of a one-parameter family of metrics along which the limit A = 0 is 
achieved, provided that those families coincide for A —> 0 and A = Apf^l 

In their paper, Charach and Malin did not put any restrictions on the 
choice of the coefficients A n , B n , C n . D n . It is rather obvious that they 
must be in some sense small: the wave part h must be a small correction 
to the background r) (e.g., p is explicitly assumed to be small). Otherwise 
their result and its interpretation is not valid. In the Green-Wald approach 
this assumption is incorporated into the framework. The introduction of A- 
dependent family of metrics allows us to define the split of the metric into 
the background g ^ and the ‘perturbation’ h in a gauge-independent manner 
(provided that transformations of coordinates do not alter the limit A —» 0). 

In addition to the reasoning outlined above, Charach and Malin sug¬ 
gested an alternative procedure to show the null fluid type high-frequency 
behavior of the scalar held source. We may average the energy-momentum 
tensor. Namely, we take the energy-momentum tensor (13) and apply the 


transformation (22). Next, we use asymptotic formulas for Bessel functions 
for £/A 1. Averaging over the phase the nontrivial components of the 

energy-momentum tensor gives 

<T K >=<T„>=A,(l + iW.„), 


^ T 


In the leading order we have 


< >=< T zz >= 


SttV 


< tr{T ) >=< T x x >=< T y y >= 0 , 


(28) 


with the remaining components equal to zero and in accord with Charach 
and Malm’s asymptotic averages: the equations (62), (63), and (64) in [6]. 
Of course, this is the averaging of components of a tensor, so it is not a 
covariant procedure. The result is invariant against a restricted class of 
coordinate transformation, but this class has not been specified in [6]. In 
order to overcome this kind of difficulties, Green and Wald averaged the 
energy-momentum tensor differently, namely, they took a weak limit of T 
(we denote it as T^ 0) = w-liniA-s-o T). Thus, in their calculation one ‘goes 


6.2 


30 Using terminology introduced in [2], there is no path dependence so far. 
we will show that the final results are also path independent. 


In Section 
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with A to zero’ along a decreasing sequence of small numbers 1/m, where 
m — j->j + l,...,+oo. Such a procedure is needed only to average the 
energy-momentum tensor T in a more gauge controlled way than in the 
Charach-Malin approach which has been demonstrated above. In the case 


under investigation the Green-Wald weak limit coincides with (28), and hence 
in the leading order < T >= T (0 h 

Charach and Matin did not discuss backreaction in their model; however, 
it seems reasonable to follow the Green-Wald approach and define backre¬ 
action in the Charach-Malin approach by the effective energy-momentum 
tensor t = T— < T >. Thus, in the Green-Wald and Charach-Malin ap¬ 
proaches the beckreaction effect for the particular choice of solution studied 
in this paper is identical. It is given by the effective energy-momentum tensor 
which is traceless and corresponds to gravitational radiation 


t = 




8vr 2 £ 


(,de + dz 2 ) = , 


(29) 


where c°) denotes the effective energy-momentum tensor in the GW ap¬ 
proach. 


6.2 Uniqueness of the effective energy-momentum tensor 

An interesting question arises: had we chosen A~, Cj differently, would the 


GW framework lead to a different effective energy-momentum tensor (29)? 


Let us assume for a moment that = g and g(Ao) = g for some large 
natural number j = 1/Ao- Physical inhomogeneities are always of finite size, 
so it seems that both conditions should define backreaction in the model 
uniquely. However, there may exist many one-parameter families g( A) that 
satisfy these conditions. The Green-Wald framework provides a mapping 3F : 
g\ —> id 0 ' as described in the Introduction, and thus backreaction for different 
one-parameter families may differ. In [2] this was called a ‘path dependence.’ 
We show below that in the model under investigation the results are ‘path 
independent.’ 


Let us consider modified solutions that differ from (23) in the choice 


A„ = , C n = VW . 3 = j; ■ 


(30) 


It follows from the definitions above that A~ 


A K / 2 and Cj ~ A^ 2 . Of course, 


other choices are also possible, but we would like to have g*- 01 = rj, g( Aq) = g 
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and since j is assumed to be large, we are interested in the asymptotic be¬ 
havior of Aj and C 3 \ thus (30) seems to be a reasonable parametrizationp] 
The solutions depend now on three more constants k, l and a. The constant 
a must be positive because otherwise A —> 0 would not correspond to the 
high-frequency limit. The leading term in the formula for / has the form 


/ « - (a 2 \ b - a + f3 2 \ K ~ a ) £ • 

7T 


(31) 


In order to recover the background g in the limit A —» 0 we must have i — a — 
k — cr = 0; thus i = k = a and there is only one independent parameter a. 
However T, = w-liuLA-s-o T and the leading terms in < T > do not depend 
on aW Thus, without loss of generality we may set cr = 1 to recover solutions 
studied in this paper. To sum up, although there may exist many one- 
parameter families g\ with the property = rj, g{ Ao) = g the Green-Wald 
framework in the model under investigation provides a unique procedure to 
calculate backreaction effects. The condition g = gA> ■= liniA^o ,9a fixes 
uniquely the asymptotics in A of m which determines the effective energy- 
momentum tensor A 0 ' 1 = T — T^°\ 33 


6.3 More general solutions 


The crucial difference between the Charach-Malin and the Green-Wald ap¬ 
proaches is illustrated by the following fact. In both approaches the effective 
energy-momentum tensor is traceless. In the Charach-Malin approach this 
is rather an ‘interesting result of calculations’ that were conducted for a par¬ 
ticular solution to Einstein’s equations. In the Green-Wald approach this is 
a consequence of their theorems (as explained in the Introduction), and thus 
a fundamental property of all effective energy-momentum tensors that may 
be calculated within their framework. 

In order to illustrate this property one may consider more general solu¬ 
tions than those studied in this paper, namely, we may take the Charach- 


31 In fact, the analysis presented in this section covers a wide class of A-dependent families 
of solutions g( A) for which the leading terms of some general functions A n { A), C n ( A), j( A) 
for small A have the form (301. 

32 Although Aj, Cj and j depend on a this dependency cancels in appropriate asymptotic 
expressions for T, < T >, T ^ by the same mechanism as in the limit 77 := lim^o gx, 
where the A approaches zero along a discrete sequence of values. 

33 Since ?? = g(°\ then T = g ^G(g^). 
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Malin solution [6] with the constant j3 0 ^ 0 and the remaining constants cor¬ 
responding to those studied in this paper. This solution contains a spatially 
homogeneous component of the scalar field and it is similar in some sense to 
the Belinskii-Khalatnikov solution. The equation (46) in [6] implies that for 
/3o 7^ 0 the energy-momentum tensor T is not traceless \tr(T) = 
and one may try to use this fact to construct a counterexample to the Green 
and Wald theorems. The equation (64) in [6] shows that for £ 1 the aver¬ 

age energy-momentum tensor < T > is traceless up to the order £ -1 . This 
is consistent with the fact that the trace of the effective energy-momentum 
tensor t = T— < T > should vanish in the highest order (terms £ -1 ) and 
remains in agreement with the Green and Wald theorems provided that we 
choose (3q ~ A. If /3o will not be proportional to A, then the trace of T may 
contribute to the trace of the effective energy-momentum tensor t. Thus, the 
Green and Wald theorems imply that there should appear a nonzero trace 
in the leading order of < T > to cancel the trace of T and make t traceless. 
Our analogy between the Green-Wald and Charach-Malin approaches sug¬ 
gests that also T ^ cannot be traceless. Indeed, calculations show that in the 
leading order tr(T) = tr(< T >) = tr(T ( '°' > ) as predicted by Green and Wald 
and tr(t) = tr(t ^) = 0 also for /3 0 7^ 0 in agreement with their theorems. 


6.4 Vacuum limit: The Isaacson approach 


The solution (|24j) corresponds for a = 0 to vacuum three-torus polarized 
Gowdy cosmologies]^] These cosmologies were already studied in a different 
coordinate system in the context of the GW framework in [40j. The results 
presented there are consistent with ours, but the direct comparison is ob¬ 
scured by nontrivial coordinate transformation (the tensor fj, a bcdef was not 
explicitly given in [ID]). The vacuum case a = 0 may be investigated within 
the Isaacson approach El- 

We find it instructive to apply the Isaacson framework to the solutions 
studied in this paper. It follows from (25) and (26) that our metric (21) 


may be written in the Isaacson form g = r) + A h, where h = h/X. We 
have h a b = 0(A°), V a hb c = 0(A -1 ) and V a Vf,/r C( i = 0(A~ 2 ), where V a is a 
covariant derivative associated with //. In the Isaacson approach the Ricci 
tensor is calculated in terms of h and its derivatives. This expression is 


34 The full metric g, which is given by (21), is a vacuum solution, but the background 


metric 77 is nonvacuum. Also 77 +h does not satisfy the vacuum Einstein’s equations exactly. 
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expanded in orders of A 


RatVaW] = d,?W + Ai&’M + A 2 i?®(A) + 0(A). 


We have A-R^(A) = 0 (A and the remaining two terms R^( A), A 2 /?^(A) 
are of the same order 0(A°). The vacuum Einstein’s equations -R a b[^(A)] = 0 
imply that i?^(A) = 0 and R^( A) + A 2 /?,^ (A) = 0; thus the average of 
the may play the role of the effective energy-momentum tensor for the 
background metric [-R^(A) = R a b(v)\- Therefore, following Isaacson | 8 | and 


Brill and Hartle [2j| (see also [12]) we dehne t ab = — ^ < R\ b ( A) >. In order 
to calculate it one may use equation (2.8) in [?], but it is more convenient 
for us to expand the final expression for R ab ( A) in powers of A, and indentify 


Rat, ( A) as R ab {jj)] the remaining terms in the order 0(A°) will give us — Af. 
We find t = t = t^ in accord with (29). Therefore, all three frameworks (by 
Isaacson, by Charach and Malin, and by Green and Wald) predict the same 
backreaction in the vacuum limit. 


7 Summary 

The one-parameter family of exact nonvacuum solutions to Einstein’s equa¬ 
tions presented in this article satisfies all assumptions of the Green-Wald 
framework [T|. The other three examples of such families presented in the 
literature so far were restricted to vacuum mm or to stiff-fluid spacetimes 
[9]. In the cosmological context, the most interesting are nonvacuum fami¬ 
lies. The novel important property of the example presented in this article 
(in comparison to the nonvacuum example presented in (9]) is the nontriv¬ 
ial behavior of the energy-momentum tensor. The matter density exhibits 
bounded nonvanishing amplitude, high-frequency oscillations. Our calcula¬ 
tions confirm the mathematical consistency of the Green-Wald approach. 
We showed that within the model studied it predicts the backreaction effects 
uniquely for a finite size of inhomogeneities. We used the family of exact 
solutions in question to compare the Green-Wald framework to the Charach- 
Malin approach, and in the vacuum limit to the Isaacson method. Although 
these methods use slightly different mathematical formalisms, all of them 
agree in their range of applications. 
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A task for the future is to find an example of a one-parameter family of 
solutions to Einstein’s equation^] for which the amplitude of fluctuations 
of matter density becomes unbounded as the background spacetime is ap¬ 
proached. 
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